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Abstract: The discoveries of the GZK cutoff with the HiRes and Auger Observatories and the discovery 
by Auger of clustering of 60 EeV ultra-high energy cosmic rays (UHECRs) towards nearby ( 75 
Mpc) AGNs along the supergalactic plane establishes the astrophysical origin of the UHECRs. The likely 
sources of the UHECRs are gamma-ray bursts and radio-loud AGNs because: (1) they are extragalactic; 
(2) they are sufficiently powerful; (3) acceleration to ultra-high energies can be achieved in their relativis- 
tic ejecta; (4) anomalous X-ray and 7-ray features can be explained by nonthermal hadron acceleration 
in relativistic blast waves; and (5) sources reside within the GZK radius. Two arguments for acceleration 
to UHE are presented, and limits on UHECR ion acceleration are set. UHECR ions are shown to be able 
to survive without photodisintegrating while passing through the AGN scattered radiation field, even if 
launched deep in the broad line region. UHECR injection throughout cosmic time fits the measured en- 
ergy spectrum of UHECRs, at least for protons. Local UHECR proton and ion interaction and energy-loss 
mean free paths are calculated using an empirical fit to the extragalactic background light (EBL) at IR and 
optical energies. Minimum intergalactic magnetic (IGM) fields ~ 10^^^ G are derived from clustering 
assuming specific source origins, e.g., Cen A, nearby AGNs, or GRBs for the super-GZK CRs seen with 
Auger. Besides distinct cosmic-ray induced 7-ray signatures that should be observed with the Gamma ray 
Large Area Space Telescope (GLAST), source and GZK neutrino detections and the arrival distribution of 
UHECR in direction and time can finally decide the sources of cosmic rays at the highest energies. 



Introduction 

A high-significance steepening in the UHECR 
spectrum at energy E = 10^^'^ eV was re- 
ported earher this year by the HiRes collabora- 
tion [1], and here at the 2007 Merida ICRC based 
on observations taken with the Auger Observa- 
tory [2]. This result confirms the prediction of 
the Greisen-Zatsepin-Kuzmin (GZK) cutoff [3, 4] 
in the UHECR spectrum due to photohadronic in- 
teractions of UHECRs with photons of the cosmic 
microwave background radiation (CMBR), and fa- 
vors astrophysical bottom-up vs. particle physics 
top-down scenarios for the UHECRs, provided that 
sources are found within the GZK radius. 

At the same time. Auger data shows [5] evidence 
for mixed composition with substantial ion con- 
tent in UHECRs with energies as high as a few x 
10^^ eV, based on studies of the depth of shower 
maxima. With hybrid fluorescence detectors and 
shower counters. Auger provides the strongest evi- 



dence yet for metals in the UHECRs, possibly with 
mean atomic mass (A) ~ 8 - 26, significantly 
different from pure proton and pure Fe composi- 
tion. This result depends on the accuracy of the 
nuclear interaction physics used to model show- 
ers, but points to the importance of nuclei in the 
UHECRs, and the meaning of this for GZK physics 
[6,7, 8]. 

The GZK radius of a proton with energy E = 
Epar = 10^°i?2o Can be estimated by not- 
ing that the product of the cross section and 
the inelasticity in a pion-producing reaction is 
K^T^a^-Pi « 70 pb [9], so that the photo- 
pion energy-loss pathlength r^T,{E) is given by 
nph{E){K^T,(74,Tr)r^T,[E) = 1, where nph{E) is 
the i?-dependent number density of photons above 
the pion-production threshold. 

There are sa 400 CMBR photons per cm'^ at the 
present epoch, and ions with Lorentz factor 7 = 
E/ArripC^ satisfying ^/{hv/rrieC^) }t 2171^^/171^ = 
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500 will interact with most of the photons 
of the CMBR. The mean dimensionless pho- 
ton energy in the CMBR is h(y)lrac(? = 
2.1Qk^TcMBRlmeC^ = 1.3 x 10"^, noting that 
TcMBR — 2.72 K at the present epoch; hence 
UHECR protons with Lorentz factor 7 ^ 4 x 10^^, 
or E' ^ 4 X 10^° eV will have a photo-pion energy- 
loss pathlength r^^{E > 4 x 10^" eV) « 12 
Mpc. By considering the number of photons above 
threshold in a blackbody distribution, the energy- 
loss mean free path of an UHECR proton with en- 
ergy E is found to be 



r07r(£'2o) = 



13.7exp[4/£;: 



20J 



[1 + 4/S20] 



Mpc . (1) 



Figure 1 shows the photopion energy-loss path- 
length from eq. (1) for cosmic-ray protons inter- 
acting with CMBR photons, in comparison with a 
more detailed calculation [10]. 

Auger results show a mixed UHECR composition 
extending from w 4.5 x 10^'' eV - 4.5 x 10^^ eV, 
with an indication of increased average mass at the 
highest energy datum [5]. The ionic content of the 
UHECR provides new information to understand 
UHECR source properties, especially if the distri- 
bution of atomic mass A can be obtained with bet- 
ter statistics. 

Figure 1 shows expansion, photohadronic, and to- 
tal energy-loss mean free paths (MFPs) for p and 
Fe on the CMBR. Calculations of photopair losses 
follow [11]. The photopion loss calculation can be 
extended to ions of atomic mass A and charge Z by 
assuming that the photopion cross-section times in- 
elasticity is at most = A-^/'^A'^ttC^tt above thresh- 
old. The photopion energy-loss rate of Fe with this 
assumption for the product of the inelasticity and 
cross section is shown in Figure 1. Because of the 
greater mass of Fe than p, far fewer photons are 
available for photopion interactions than with pro- 
tons carrying the same energy, so that the corre- 
sponding GZK photopion radius is much larger 

Ions are in addition subject to losses due to pho- 
tonuclear interactions that can break up nuclei (for 
example, via the giant dipole resonance). The pho- 
todisintegration loss length for Fe in the CMBR 
at z ^ 1, calculated from model fits to photo- 
hadronic cross section data [12, 13], is plotted in 
Figure 1. 
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Figure 1 : Energy-loss mean-free path for UHECR 
protons (black) and Fe ions (gray) in the CMBR 
at 2: = 0, as a function of total particle en- 
ergy. The short-dashed black curve give the ana- 
lytical approximation for the UHECR proton pho- 
topion energy-loss pathlength, eq. (1), in compar- 
ison with the numerical results, shown by open 
circles [10]. Numerical results for the photopair 
energy-loss pathlength of UHECR protons and Fe 
are shown by the dotted curves. The gray dot- 
dashed curve gives the photodisintegration energy- 
loss curve of Fe. Expansion losses are shown by 
the long-dashed line, with c/Hq = 4170 Mpc. 



UHECR Fe is seen to have a comparable GZK- 
type cutoff, but here at w 2 x 10^" eV, with pho- 
topair losses playing a minor role. 

A realistic calculation of UHECR ion spectra in 
the evolving background radiation fields must fol- 
low a reaction pathway using Monte Carlo tech- 
niques, with the actual effective energy-loss path- 
length taking into account the EBL intensity con- 
tributed by stars and black holes, then sometimes 
reprocessed through dust and gas. The exact form 
of the EBL intensity between « 1 and 100 fi is 
poorly known, but can be constrained by empirical 
galaxy SEDs and 7-ray observations. 

From the curves shown in Figure 1, we see that 
for either UHECR protons or Fe ions, the sources 
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of E ^ 10^° eV cosmic rays must be found within 
a few hundred Mpc, and the 3 x 10^° eV cosmic 
ray detected by the Fly's Eye [14] must have origi- 
nated from a source within a few tens of Mpc. Both 
gamma-ray bursts (GRBs) and radio-loud active 
galactic nuclei (AGNs') satisfy this requirement. 
Besides the availability of extragalactic sources 
within tens of Mpc, several other conditions must 
also be met, including 

1 . sufficient power to make the UHECRs at the 
observed intensity; 

2. a plausible mechanism to accelerate parti- 
cles to ^ 10^° eV energies; and 

3. a model that can reproduce the observed 
UHECR spectrum, consistent with astro- 
nomical properties of the sources. 

The Auger clustering observations [15] suggest 
new calculations to determine the intergalactic 
medium (IGM) magnetic field. If the 2 (out of 27 
total) UHECRs with arrival directions within 3° of 
Cen A are hypothesized to originate from Cen A, 
then lower values of the mean IGM field can be 
derived, ^ 10~^° G, between us and Cen A, and 
^ 10~^^ G for sources at 75 Mpc. The cluster- 
ing observations and the GZK length are tightly 
coupled concepts [16], as is established with the 
Auger data. Here we give new proton and ion en- 
ergy loss calculations after revisiting the problem 
of theEBL. 

The Auger discovery [15] that the arrival directions 
of > 56 EeV are correlated with nearby, <, 75 
Mpc AGNs, in particular, 2 cosmic rays within 
3° of Centaurus A, does not mean that AGNs are 
the sources of UHECRs, any more than O and B 
stars correlated with galactic 7-ray sources mean 
that cosmic rays are accelerated by high-mass stars 
[17]. Long duration GRBs trace sites of active star 
formation, and may be correlated to AGN activ- 
ity. AGNs trace the local matter distribution due 
to structure formation, as do radio galaxies and 
GRBs, or for that matter clusters of galaxies and 
radio-quiet Sy AGN. Here we make the case for 
GRB and radio-loud AGN/blazars as the sources 
of the UHECRs. UHECR ion acceleration in clus- 
ters of galaxies could also make a contribution at 
^ 10^^ eV, or from dim quasars [18], but this anal- 



ysis forces us to reject UHECR origin from radio- 
quiet AGNs. 

Effects of UHECR hadron acceleration can be 
varied, including anomalous 7-ray emissions 
and characteristic behaviors of GRB X-ray light 
curves. Detection of high-energy neutrinos from 
hadronic interactions at the sources will provide 
the most definitive evidence for the presence of 
UHECRs [19]. Detection of cosmogenic GZK 
neutrinos from the decay of pions formed in GZK 
interactions with CMBR photons gives a calori- 
metric measure of UHECR power throughout cos- 
mic time, with a characteristic spectrum imprinted 
by production and propagation effects. ANITA and 
successor Askaryan-effect detectors should soon 
be able to constrain a long-duration GRB model 
of UHECRs, if they were protons. Metals in the 
composition of UHECRs change the spectral pre- 
dictions, as well as the GZK neutrino predictions 
[20, 21]. The result of UHECR ions vs. protons is 
mainly to enhance neutron /3-decay neutrinos but 
not significantly change the ^ EeV neutrino flux. 

Joint analysis of UHECR spectral, composition, 
and directional information, now given most ac- 
curately with the hybrid Auger, with neutrino and 
7-ray data from GLAST, AGILE, and the ground- 
based air and water Cherenkov detectors, should 
lead to a definitive solution to the problem of 
UHECR origin. A black hole origin is explored 
here. Crucial to progress are the advances in de- 
tector capabilities and instrumentation, multiwave- 
length observations and multi-disciplinary analy- 
ses. 



Extragalactic Origin and Source Ener- 
getics 

The Larmor radius of an UHECR ion [22] is 

QB ~ (Z/10)S(/xG) P 
(E/6X lO^^eV) 

where B = 10~'^^i3_ii G is the mean magnetic 
field in which the ion propagates. A significant 



1. Blazars are radio-loud AGNs viewed along the jet 
axis; here the two terms are used interchangeably. 
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Table 1: UHECR Source Emissivity 
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"From Auger data [2]. 



anisotropy in the arrival directions of UHECRs 
with energies exceeding k, lO^^ eV might be ex- 
pected from Galactic sources along the Galactic 
plane or towards the Galactic Center, even if they 
were composed primarily of Fe {Z = 26), because 
any likely source class (pulsars, supernovae, mi- 
croquasars, high mass stars) would be confined to 
the thick disk of the Galaxy. HiRes found no evi- 
dence for small-scale or large-scale clustering [23]. 
The Auger discovery [15] of UHECR clustering 
along the supergalactic plane, anticipated in anal- 
ysis of AGASA, Haverah Park, and Yakutsk data 
[16], immediately rules out a galactic origin of the 
UHECRs. 

We now calculate the emissivity £ required to 
power UHECRs, assuming (incorrectly) that UHE- 
CRs are protons that suffer expansion cooling and 
photopion and photopair losses on the CMBR. The 
Auger observations [2] of the spectrum of UHE- 
CRs, written in the form J24 = E'^J/(W^ eV^ 
vcT^ s~i sr~i), where J is the UHECR number in- 
tensity, implies that the energy density of UHECRs 
with energy £"20 is Uuhecr = 6.7 x 10"^^ J24/.E20 
ergs cm~'^. Using the energy-loss pathlengths 
shown in Fig. 1 gives^ £ = Uuhecr/ttot = 

CUuhecr/rtot, SO 



ergs 



6.0 X 10^5 J24 
rtot (Mpc) W 



(3) 



20 



The sources of UHECRs with energy « 102°i;2o 
eV must provide an emissivity £44, in units of 10'*'' 
ergs Mpc~^ yr~i, given by the values shown in 
Table 1. 

The required emissivity, from Table 1 using the 
Auger intensity, is « few x 10^^ ergs Mpc~^ 
yr-i for E > eV, and « IC^^ ergs Mpc'^ 
yr-i for E > 10^^ eV. Between lO^s eV and lO^o 
eV, the increased energy density in UHECRs at 



the lower energies of this range is balanced by the 
larger energy-loss length, thereby making the in- 
jection emissivity roughly constant with energy. At 
E ^ 10^^ eV, the emissivity increases roughly pro- 
portional to energy because the energy-loss path- 
length is approximately equal to the Hubble radius 
at IQi^ eV, and £^,7 oc 1/i? at these energies. 

Classical long duration GRBs have a volume- and 
time-averaged emissivity in the X-ray/soft 7 ray 
("X/7") band comparable to this value, a coin- 
cidence first noted by Vietri [24] and Waxman 
[25]. We can reproduce this estimate by noting 
that the average > 20 keV fluence of BATSE (the 
Burst and Transient Source Experiment (BATSE) 
on the Compton Observatory) GRBs is w 10~^ 
ergs cm~^ [26], and that there are about 500 long- 
duration BATSE GRBs over the full sky per year 
[27], giving an average GRB flux of ~ 10~^ ergs 
cm^^ yr^i.^ BATSE GRBs are, on average, at 
redshift z « 1 (Hubble radius Rn ^ 4200 Mpc), 
so that the ^ 20 keV emissivity of long duration 
GRBs is 



^GRB 



AnR^ X 10 ^ ergs cm ^ yr ^ 
*^ 47ri?3/3 

7 X 10"*^ ergs Mpc"^ 



yr ^ , (4) 

in rough agreement with the UHECR emissivity re- 
quirements shown in Table 1. 

Two assumptions (at least) underlay this estimate: 
One is that an average emissivity can apply to 
the pecuhar emissivity of the local ^100 Mpc 
sphere that we inhabit. Greater GRB activity at 
z ^ 1 compared to the present epoch means that 
there would be fewer sources within the GZK ra- 
dius [28]. On the other hand, additional classes of 
GRB sources, in particular, the X-ray flashes or the 
low luminosity GRBs [29], can provide a substan- 
tial addition to the emissivity. Wang et al. (2007) 
[30] estimate a local emissivity £ w 250 x 10^^ 
ergs Mpc"'^ yr~i in the kinetic output of nearby 
low-luminosity GRB hypemovae (which exceeds 



2. This is the proton GZK radius. Ion energy-loss 
pathlengths are calculated later in the paper. 

3. Gonzalez (private communication, 2003), calcu- 
lates a total fluence per year of 0.63 x 10^^ ergs cm^^ 
for 1293 GRBs in the 4th BATSE catalog, including 
GRBs of both long- and short-duration, considering 666 
GRBs/yr full sky, implying a bolometric average energy 
density ^ 10"^° ergs cm~^ of GRB light. 
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Figure 2: Intensity of extragalactic background 
light from the microwave through the 7-rays, with 
components as labeled. 

the emisivity in observed 7 rays). A further as- 
sumption is that there is comparable energy in- 
jected in UHECRs as detected in electromagnetic 
radiation. This could well be wrong. A large, 
« 30 - 100, baryon loading is required if long du- 
ration BATSE/Beppo-SAX/GBM-type GRBs are 
the sources of UHECRs [31]. 

A similar emissivity estimate for blazar AGNs can 
be made on the basis of the unresolved/diffuse 
extragalactic 7 radiation (DEGR) observed with 
the Energetic Gamma-Ray Experiment Telescope 
(EGRET) on the Compton Observatory [32, 33] 
(see Fig. 2). Between sa 100 MeV and « 100 
GeV, the "vFy" intensity e/e « 1.5 keV cm^^ s^^ 
sr^^, so that the diffuse 7-ray energy flux received 
here at Earth is « 0.5 erg cm^^ yr^^. Because 
blazar AGNs comprise the largest number of iden- 
tified EGRET sources, they undoubtedly make up a 
large fraction of this radiation, with estimates rang- 
ing from w 20% to nearly 100%. 

It is important here to distinguish between the 
nearly lineless BL Lac objects and the flat spec- 
trum radio quasars (FSRQs) with strong broad, op- 
tical emission lines. The BL Lac objects detected 



at 7-ray energies are at mean redshift (z) w 0.3, 
whereas the FSRQs detected with EGRET are at 
(z) w 1. The BL Lac and FSRQ contributions to 
the DEGR are found [34] from model fits to the 
redshift and peak flux 7-ray data to contribute at 
the ^ 5% and ~ 15% of the total diffuse 7-ray 
intensity, respectively. 

Following the reasoning leading to eq. 4, the 
mean emissivity in 7 rays is £bl La 



lO''^^ ergs Mpc yr ^ for BL Lac objects, and 



Sfsrq ~ 6 X lO'*" ergs Mpc"^ yr"! for FSRQs. 
On the basis of energetics arguments, both FSRQs 
and BL Lacs appear to have sufficient energy to 
power the UHECRs. The more powerful FSRQs 
are, however, rare in our local vicinity, so that av- 
erage emissivity becomes a mistaken concept. The 
nearest high luminosity FR II radio galaxies asso- 
ciated with FSRQs are at z « 0.1; the FR II radio 
galaxy Cygnus A has z k, 0.056, so is w 250 Mpc 
distant, and is far outside the GZK radius for a lO^'' 
eV proton. 

By comparison, many low-luminosity FR I radio 
galaxies associated with BL Lac objects are found 
nearby. For example, Centaurus A and M87 are 
w 4 and ~ 16 Mpc distant, respectively. These FR 
Is have lower jet power than the powerful FR lis. 
This makes it more difficult to accelerate protons 
to ultra-high energies, though not heavy ions, as 
we now show. 



Acceleration to Ultra-High Energies 

Acceleration in relativistic blast waves can pro- 
ceed through a number of mechanisms. First- 
order shock Fermi acceleration, analogous to the 
standard particle acceleration mechanism in super- 
nova remnant shocks, is the obvious process [35]. 
But acceleration of particles to ultra-high energies 
through first-order processes at a relativistic exter- 
nal shock encounters kinematic difficulties to reach 
ultra-high energies [36, 37]. When particles diffuse 
upstream from the shock, the shock overtakes the 
particles before they can change their dkection by 
more than an angle « 1 /F, where F is the Lorentz 
factor of the blast wave. This prevents the parti- 
cle from increasing its energy by more than a fac- 
tor w 2 in each cycle following the first particle 
capture by the blast wave. Even so, acceleration 
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to UHE through first-order relativistic shock accel- 
eration can take place if the surrounding medium 
is strongly magnetized, as could be the case in 
the winds of a massive Wolf-Rayet progenitor. If 
the colliding shells of blazars and GRBs are suf- 
ficiently magnetized, first-order processes can also 
apparently accelerate protons to UHEs [25]. 

Here we present two arguments showing that 
UHEs can be achieved in the relativistic shocks 
of GRBs and blazars. The first considers stochas- 
tic acceleration in the external shock of a GRB 
or blazar, and the second considers general source 
power requirements [38] for particle acceleration 
in GRBs, blazars, or other candidate UHECR 
sources. 

A. Stochastic Acceleration in tlie Blast 
Wave Shell 

Second-order processes in an external shock of a 
blast wave can be shown to permit hadronic accel- 
eration to ultra-high energies, in comparison with 
radiative and expansion losses and escape [39]. 
Using parameters appropriate to a GRB, consider 
a blast wave with apparent isotropic energy release 
Eq = 10^'*£'54 ergs, (initial) coasting Lorentz fac- 
tor Fq = SOOFsoo, and external medium density 
ng = 100n2 cm~'^. The comoving blast-wave vol- 
ume for a spherically symmetric explosion at dis- 
tance X from the center of the explosion is 



deceleration radius [40] 



(5) 



where primes refer to the comoving frame, and the 
shell width A' = x/UT (the factor 1/12T is the 
product of the geometrical factor 1/3 and the fac- 
tor l/4r from the compression of the material by 
the shock, in accord with the conservation laws of 
relativistic hydrodynamics). 
A necessary condition to accelerate to some en- 
ergy E'^g^^ is that the particle Larmor radius is 
less than the size scale of the system [22], that is, 
r'l^ ^ Xd/12r . E„iax in the stationary explosion 
frame is then given by 



E' 

max 

ZeB' 



Er, 



VZeB' 



< A'. 



(6) 



The largest particle energy is reached at the decel- 
eration radius x = Xd when F = Fq, where the 



Xd = (- 



U/3 



2.6 X 10i6(_f:51_)V3 cm . (7) 



^47rrgmpno ' 

^300^2' 

Hence E,nax = ZeB'xd/12 . 
The mean magnetic field B' in the GRB blast wave 
is assigned in terms of a magnetic field parame- 
ter €B that gives the magnetic field energy density 
in terms of the energy density of the downstream 
shocked fluid, so 



B' = (327rnoeBmpC^)^/Vr(r - 1) = 
0.47ii^F ^ 1200V^r3ooG. 



(8) 



Thus 



Ema. = 8 X 10^°Znl/^e'J^TlilEli^ eV (9) 

[41, 39], so that external shocks of GRBs can ac- 
celerate particles to ultra-high and, indeed, super- 
GZK energies. The highest energies are more eas- 
ily reached for ions of large Z, which helps to relax 
hmits on density and e^. 

For values appropriate to blazar AGNs, acceler- 
ation to UHEs in the external shock of a blazar 
plasma jet seems possible from eq. 9, at least for 
FSRQs. For this blazar class, superluminal motion 
observations [42] and constraints from the require- 
ment that the emission region be optically thin to 
77 pair production attenuation [43] imply F ^ 10 
- 30. The brightest blazar flares observed with 
EGRET reach « 10^^ ergs s^^ and last for <, 1 
day [44], suggesting again that the total isotropic 
energy release £54 « 1. For F w 30, the decel- 
eration time td ~ {1 + z)xd/T^c is shorter than 
a day, provided that the blazar blast wave passes 
through a medium of density w 10^ cm~^. With 
these numbers, UHECRs can be accelerated in the 
external shocks of FSRQs. Much improved data 
for this estimate will be provided with GLAST. 

The case for UHECR acceleration in BL Lac ob- 
jects is not as favorable. Although spectral mod- 
eling suggests that F could be as large as 50 [45], 
values of F inferred from superluminal motion ob- 
servations rarely exceed 5-10 [46]. The mean 
distances, energy fluxes, and flare durations are 
generally smaUer for BL Lacs than quasars, and 
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the total energy in BL Lac 7-ray flares is smaller 
than for FSRQ flares. The deconvolved apparent 
7-ray luminosity from PKS 2155-304 flares dur- 
ing the 2006 July flaring state was w 3 x 10^^ erg 
s~^ [47J, and lasts between rs 10^ and 10^ s, giv- 
ing a total apparent energy releases fa 3 x lO"*^ 
- 3 X lO'^* ergs. The surrounding medium den- 
sity is lower in BL Lacs than in FSRQs, given that 
blazars have smaller optical emission-line equiva- 
lent widths than FSRQs, and therefore probably a 
smaller column density of broad Une region ma- 
terial. All these factors make it harder for low-Z 
ions to reach UHEs. So, if FR I radio galaxies are 
the sources of UHECRs, then they would have a 
suppressed proton and low-Z ion content due to 
acceleration limitations. 

Why, however, consider stochastic acceleration 
processes in the blast waves formed by the exter- 
nal shocks in blazars and GRBs, given that first- 
order Fermi acceleration through colliding shells, 
which can operate in both blazars and GRBs, could 
accelerate the cosmic rays to high energies? The 
problem here is that the material which forms 
the outflowing plasma winds is processed through 
extreme environments around supermassive black 
holes of blazar AGNs and newly formed black 
holes in GRBs. Whether ions can survive, as re- 
quired by the Auger data [5], is an open ques- 
tion. If they do survive, then Wang, Razzaque, 
& Meszaros [48] have shown that acceleration of 
UHECR ions to 10^° without significant losses 
to photodisintegration in GRB internal shocks, ex- 
ternal shocks, and hypemovae is possible (see also 
[21]). 

In a colliding shell scenario, the ejecta originates 
from the vicinity of a black hole. For the fire- 
balls that power a GRB, heavier nuclei are broken 
down into protons and neutrons, deuterium and a 
particles [49]. Nuclear breakdown reactions, ei- 
ther through direct particle spallation processes or 
through photodisintegration from the intense radi- 
ation field in the vicinity of a black hole, make it 
questionable if baryonic material ejected from the 
central engines of black holes in GRBs or blazars 
has any metal {A > 4) content. 

In a coUiding shell scenario, successive waves of 
ejecta are supposed to form the flares in blazars 
and GRBs, and if the hadronic content in the ejecta 
were composed entirely of protons and Ught nuclei. 



then such an acceleration scenario could not ac- 
count for the Auger results. Wang, etal. [48] argue, 
however, that mixing instabilities could develop 
to seed the relativistic ejecta with heavier ions as 
the blast wave passes through the the stellar enve- 
lope surrounding a GRB. In this case, the amount 
of baryon contaminant would have to be carefully 
regulated without heavily loading and quenching 
the fireballs, unless dirty fireball bursts [50] and 
quenched, or choked [51] fireballs also occur 

A scenario involving an external shock would 
permit the capture of ions from the surrounding 
medium. This medium could be highly enriched 
from the presence of circumnuclear starbursts sur- 
rounding the supermassive black hole in a blazar 
AGN or OB associations in which a GRB might re- 
side, and so have a considerable ion content. Even 
if the material had sub-Solar metallicity, heavier 
ions would be preferentially accelerated due to 
their larger charge Z (see eq. 9), producing a mixed 
composition in the UHECRs. 

B. Power Requirements for Electrodynamic 
Acceleration 

We apply Waxman's argument [38], giving mini- 
mum source power to accelerate 10^° £'20 eV pro- 
tons, to ions. In a region of size R and mag- 
netic field B, electromagnetic forces can acceler- 
ate a particle to a maximum energy of Emax > 



E„ 



ZePBR. The available time in the co- 



moving frame is shortened by bulk Lorentz fac- 
tor r, so that the effective size for acceleration is 

w R/T, and BR > TEpar/Ze/S. The required 
power of the magnetized flow is 

r2 cT^ E"^ 

L^2xAnR\xl-^PciBRr^^^ 



or 



^ 3 X 10^5 r2 2 , 



5 X 10*2 r2 



(Z/26)2 f3 



Eio ergs s ^ , 



(10) 



including a factor of 2 for the plasma jet kinetic 
power. If the nonthermal luminosity is a good mea- 
sure of jet power, then we can decide whether dif- 
ferent source classes are good candidate UHECR 
sources. 
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For GRBs, T w 300, and Lgrb » 
1050(r/300)2£;2^/Z2 ergs s'K Apparent 
isotropic X/^ powers of long-duration GRBs 
are regularly measured in excess of 10^° ergs s~^ 
|52|, so long-duration GRBs are a viable candidate 
for UHECR acceleration. For the low-luminosity 
GRBs, which may only reach « 10''* - 10^° ergs 
s~^ [53], higher Z ions can stiU be accelerated 
to super-GZK energies if T remains large. It will 
be interesting to compare radiative powers with 
this minimum power using allowed values of F 
determined from GLAST data through 77 opacity 
arguments, for different classes of GRBs. 

For blazars with F = 10, Lhiazars ^ 3 x 
10*^(r/10)2£:|o/Z2 ergs s'^. FSRQ blazar 7- 
ray flares brighter than ~ 10^^ ergs s~^ were fre- 
quently detected with EGRET [54]. Present data 
does not exclude FSRQ blazars from being the 
sources of UHECRs, especially if the accelerated 
UHECRs are primarily heavy ions. Comparisons 
between measured 7-ray luminosity and minimum 
power requirements using values of F obtained 
from 77 arguments, both of which can be deter- 
mined from GLAST LAT data, can decide whether 
FSRQ blazars can accelerate UHECR protons and 
ions. 

Eq. (10) shows that it is more difficult to acceler- 
ate UHECRs, especially UHECR protons, in the 
lower luminosity X-ray/TeV blazars with apparent 
7-ray powers ;S 10*^ ergs s~^. If TeV observa- 
tions with VERITAS, HESS, or MAGIC were to 
require, either from spectral modeUng or 77 argu- 
ments, tiiat F > 50 in sources like Mrk 421, Mrk 
501, or PKS 2155-304, acceleration even of Fe to 
super-GZK eneriges might be problematic in the 
BL Lac sources. 

It is also interesting to apply eq. (10) to parameters 
of merging clusters of galaxies, which have also 
been studied as a potential source of UHECRs [55, 
56]. 

The gain in kinetic energy when the minor cluster 
of mass M2 = lO^^Mi4M0, treated as a test parti- 
cle in the mass distribution of the dominant cluster 
of mass Ml = IO^^Mi^Mq, falls from radius ri 
to radius r2(< ri), is M2i;|/2 = GMiM2(r^^ - 
rf ^) (e.g., [57]). Thus v ^ ^/2GMJr^ when 



r2 -C ri, and 



V2 = I3cic « 3000 



Mis 



r2(Mpc) 



km s 



(11) 



so Pel — 10~^. If r2 is scaled to a typical core 
radius of the dominant cluster, ^ 0.25 Mpc, the 
greater power output during this merging episode 
occurs during a timescale shorter by a factor r'^/^. 
This can be shown by noting that the characteris- 
tic merger time t is determined by the accelera- 
tion a = GMi/r\ at the outer radius. Because 
n w a2p/2, 



/ 2rl 
GMi 



660 



3/2 

Mi5 



Myr , (12) 



where ri = rjvipc Mpc. The available energy in the 
collision is 



GM1M2 8 X 10^3 



n 



^Mpc 



Mi5 Mi4 ergs . 



(13) 

The ratio of eqs. (13) and (12) gives the maximum 
power available from the gravitational energy of 
the merging clusters, namely 



Lci = ^=4:X 10"^ 
t 



Ml 
r 



5/2 

Mpc 



Mi4 ergss-i, (14) 



When rMpc ~ 0.25, corresponding to typical core 
radii of galaxy clusters like Coma, the maximum 
merger power is w 30 larger and the timescale, eq. 
(12), is a factor « 8 less. The luminosity require- 
ment to accelerate UHECRs with energy is, 
for parameters typical of merging cluster, 

V ^ 3 X 10*«i;|o/[^'(/3d/10-')] ergss-1 , 



from eq. (10). According to this criterion, it is not 
difficult to accelerate UHECR protons in galaxy 
cluster environments. This estimate does not con- 
sider actual timescales [58] of acceleration, which 
result in maximum proton energies ^ 10^^ eV for 
p from nonrelativistic shocks [59]. Acceleration of 
p to < 10^9 eV and Fe to < 10^° eV is possible 
[56] in large Mach number [60] cluster accretion 
shocks. 

The time for the merger during maximum power 
output, from eq. (12), corresponds to a distance 
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f = ci ^ 25(rMpc/0.25 Mpc)^''/VAfi5 Mpc, 
and UHECR Fe might go through some significant 
photo-erosion in the CMBR and EBL, so Fe would 
have difficulty surviving to ^ 10^° eV if it were 
accelerated by merger shocks in merging clusters 
of galaxies (cf. [58] for cluster accretion shocks). 

X-ray and 7-ray Signatures of UHECR 
Acceleration 

Indirect evidence for UHECR acceleration is given 
by analysis of spectra and light curves of GRBs and 
blazars. In the relativistic blastwave framework, 
the hard X/j radiation during the prompt phase of 
a GRB is primarily nonthermal synchrotron radia- 
tion, including perhaps some thermal photospheric 
emission [61]. A synchrotron self-Compton (SSC) 
component would accompany the nonthermal syn- 
chrotron emission, and display a correlated behav- 
ior. 

The delayed hard emission tail in GRB 940217 
[62], lasting for over 90 minutes after the start of 
the GRB, could be a signature of a long lasting, 
hadronic acceleration process. The slower decay of 
the hadronic emission component compared to the 
leptonic component, as expected from the standard 
blast wave model [63], might explain the delayed 
emission. Nevertheless, a leptonic model with 
an SSC component appearing in the GeV regime 
when the synchrotron component has decayed to 
optical/UV energies could also account for the de- 
layed behavior of GRB 940217, or the superbowl 
burst, GRB 930131 [64]. 

A more difficult case to explain in terms of non- 
thermal lepton radiations is GRB 941017, ob- 
served with the BATSE Large Area Detector 
(LAD) and the EGRET Total Absorption Shower 
Counter (TASC). Gonzalez et al. [65] reported the 
detection of an anomalous MeV emission compo- 
nent in the spectrum of this burst that decays more 
slowly than the prompt emission detected with the 
LAD between w 50 keV and 1 MeV range. The 
multi-MeV component lasts for > 200 seconds, 
and is detected both with the BATSE LAD near 
1 MeV and with the EGRET TASC between w 1 
and 200 MeV. The spectrum is very hard, with a 
photon number flux ^{e^) oc e~^, where e-y is the 
observed photon energy. 



This component was not predicted nor is easily 
explained within the standard leptonic model for 
GRB blast waves. It has been suggested that this 
emission component could be related to Comp- 
tonization of reverse-shock photons by the forward 
shock electrons, including self-absorbed reverse- 
shock optical synchrotron radiation [66]. Ex- 
tremely large apparent isotropic energies are how- 
ever required. 

This component could be a consequence of the 
acceleration of ultrarelativistic hadrons at the rel- 
ativistic shocks of GRBs [67]. A pair-photon 
cascade initiated by photohadronic processes be- 
tween high-energy hadrons accelerated in the GRB 
blast wave and the internal synchrotron radiation 
field produces an emission component that ap- 
pears during the prompt phase, as shown in Fig. 
3, but delayed due to the time required for ac- 
celeration. Photomeson interactions in the rel- 
ativistic blast wave would simultaneously make 
a beam of UHE neutrons and neutrinos, as pro- 
posed for blazar jets [68]. Subsequent photopion 
production of these neutrons with photons out- 
side the blast wave will produce a directed hyper- 
relativistic electron-positron beam in the process 
of charged pion decay and the conversion of high- 
energy photons formed in 7r° decay. These ener- 
getic leptons produce a synchrotron spectrum in 
the radiation reaction-limited regime extending to 
^ GeV energies, with properties in the 1 - 200 
MeV range similar to that measured from GRB 
941017. GRBs displaying anomalous 7-ray com- 
ponents are most likely to be detected as sources of 
high-energy neutrinos with IceCube. 

If UHECRs are accelerated by GRB blast waves, 
then blast- wave dynamics will be affected by the 
loss of internal energy when the UHECRs escape. 
This effect [69] could explain the rapid X-ray de- 
clines in the Swift GRB light curves [70]. Protons 
undergoing photopion interactions with photons at 
the peak photon frequency or peak dimension- 
less energy epk = hUpk/meC^ = 2Te'p^/{l + z) of 
the vF,^ spectrum have energy, as measured by an 
observer in the explosion frame, 

Ep, = mpc%, ^ (1 + ,),^, ^V. (15) 

The comoving time required for a proton with en- 
ergy Epk to lose a significant fraction of its energy 




Figure 3: Atoyan's calculation of photon energy 
flux from an electromagnetic cascade initiated by 
photopion secondaries in a model GRB, with pa- 
rameters given in Ref. [67]. Five generations of 
Compton (heavy curves) and synchrotron (light 
curves) are shown. The first through fifth gener- 
ations are given by solid, dashed, dot-dashed, dot- 
triple-dashed, and dotted curves, respectively. The 
total cascade radiation spectrum is given by the up- 
per bold dotted curve. 

through photohadronic processes is 

^ 2 X 10« ^MV52)ii±f)fl^ 3 , (16) 

"28/- 6 

where x = lO^^xie cm and /c^^ — 10^^/-6 ergs 
cm~^ s^^ is the lyF^ flux measured at epk', the re- 
lation between Epk and epk is given by eq. (15). 

The dependence of the terms x{t), fe ^{t), r{t), 
and Cpkit) on observer time in eq. (16) can be ana- 
lytically expressed for the external shock model in 
terms of the GRB blast wave properties Eq, Tq, en- 
vironmental parameters, e.g., no, and microphys- 
ical blast wave parameters es and Ce [69]. This 
can also be done for other important timescales, 
for example, the (available) comoving time t'^.^^ 
since the start of the GRB explosion, the comov- 
ing acceleration time t'^^^ — C,ac.cmpC^l' / cBc, 
written as a factor C,ac.c ^ 1 times the Larmor 
timescale [71] (Cacc = 10 in the Figure 3), the es- 
cape timescale t'^^^ in the Bohm diffusion approx- 
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Figure 4: Rates and inverse timescales as a func- 
tion of observer time for 10'^" eV cosmic ray pro- 
tons as measured by a stationary external observer. 
The figure uses parameters given in the text. 

imation, and the proton synchrotron energy loss 
timescale t'^y^. 

Fig. 4 shows the rates (or the inverse of the 
timescales) for 10^" eV protons in the case of an 
adiabatic blast wave that decelerates in a uniform 
surrounding medium. The parameters are 

2 = 1, To = 150, E^i = 10, no = 1000 cm"^ 

ee — 0.1, and = 0.3 . 

For these parameters, it takes a few hundred sec- 
onds to accelerate protons to energies sa 10^" eV, 
at which time photohadronic losses and escape 
start to be important. Photohadronic losses inject 
electrons and photons into the GRB blast wave. 
The electromagnetic cascade emission, in addi- 
tion to hyperrelativistic electron synchrotron radi- 
ation from neutron escape followed by subsequent 
photohadronic interactions [67], makes a delayed 
anomalous 7-ray emission component as observed 
in GRBs 940217 and 941017 [62, 65]. Ultra-high 
energy neutrino secondaries are produced as by- 
products of the hadronically-induced cascade. The 
ultra-high energy neutrons and escaping protons, 
accompanied by escaping UHECR ions, form the 
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UHECRs with energies > 10^° eV. Detection of 
high-energy neutrinos from GRBs would confirm 
the importance of hadronic processes in GRB blast 
waves. 

The GRB blast wave quickly loses internal energy 
due to the photohadronic processes and particle es- 
cape. The blast wave will then decelerate, produc- 
ing a rapidly decaying X-ray flux. I [69] argue that 
the decaying fluxes in Swift GRBs are signatures 
of UHECR acceleration by GRBs. (See [72] for 
other explanations of the Swift data.) If this sce- 
nario is correct, GLAST will detect anomalous 7- 
ray components preferentially in those GRBs that 
undergo rapid X-ray declines in their X-ray light 
curves. 

Anomalous 7-ray signatures have also been de- 
tected in the spectra of blazar AGNs, for exam- 
ple, the orphan 7-ray flare observed in the TeV 
blazar lES 1959H-650 [73]. This is a case where the 
correlated variability between the synchrotron X- 
rays and SSC 7 rays expected in the standard syn- 
chrotron/SSC TeV blazar model is not observed. 
Orphan X-ray flares from an hadronic emission 
component could be produced by cosmic ray pro- 
tons with Lorentz factors 7 w 10^ - 10'' under- 
going photohadronic interactions with reflected X- 
ray target photons [74]. The implied neutrino sig- 
nature is unfortunately too weak to be detected 
with IceCube [75] or a Northern Hemisphere km- 
scale neutrino telescope. In order for a TeV 
blazar to have efficient photohadronic interactions 
of UHECR protons with ambient jet synchrotron 
radiation, the blazar must be optically thick to 77 
pair production [76, 77], so we would only expect 
PeV neutrinos from TeV blazars during times of 
low-state TeV 7 ray flux. Important for this search 
is multiwavelength GRB and blazar capability. 

In the case of FSRQ blazars, there is as yet 
no strong evidence for anomalous 7-ray emission 
components that could be associated with UHECR 
acceleration. This lack of evidence should not be 
considered definitive, for a number of reasons: (1) 
The leptonic models for FSRQs are more compli- 
cated than for TeV blazars or GRBs, and involve 
a variety of external radiation fields and, conse- 
quently, more parameters. Back-scattered radia- 
tion in structured blazar AGN environments makes 
another radiation feature. This makes it more diffi- 
cult to ascribe an emission component to a non- 
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Figure 5: Different star formation rate (SFR) his- 
tories, as described in the text. 



leptonic origin. (2) Except in a few cases, e.g., 
[78, 79], the 7-ray data from FSRQs taken dur- 
ing the EGRET era didn't have extensive multi- 
wavelength monitoring. (3) The sensitivity of the 
EGRET telescope permitted spectral fits with a sin- 
gle power law integrated over several days in all 
but a few cases. This situation will change dra- 
matically with the launch of GLAST in 2008, and 
is already changing with the development of low- 
energy threshold air-Cherenkov 7-ray detectors. A 
significant advance in this direction was reported 
at this conference by the MAGIC collaboration 
[80], namely the 6 and 5 a detections of 3C 279 
(z = 0.538) in the respective bands ^ 80 - 220 
GeV and ^ 220 GeV This opens the question of 
what other nearby FSRQs will be detected with 
MAGIC, and what this means for the intensity of 
the EBL. 

UHECR Protons from GRBs 

GRBs are argued to be the sources of the UHE- 
CRs for over a decade [24, 25, 81]. Well-defined 
calculations based on particle physics and GRB as- 
tronomy for UHECR proton propagation in the ex- 
panding universe, treating expansion, and photo- 
pion and photopair interactions with CMBR pho- 
tons only (see Figure 1), give results that can be 
compared with data and used to benchmark more 
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Figure 6: UHECR data from monocular HiRes 
(open crosses and circles), Auger (filled circles), 
and AGASA (filled diamonds) data, compared 
with GRB predictions for UHECR protons with 
different SFRs shown in Fig. 5. Protons are in- 
jected from 10^'' eV to an exponential cutoff en- 
ergy E„iax = 10^" eV with local emissivity = 
60 X 10*'' ergs Mpc"^ yr"^ 



detailed calculations involving nuclei and various 
source classes. 

Restricted to the long-duration class, of which we 
have the most knowledge, the act of making a cal- 
culation of the UHECR spectrum from GRBs re- 
quires (i) an injection spectrum, which we take 
to be a power-law with number injection index 
s = 2.2 between lO'^'' eV and an exponential cut- 
off energy, usually taken to be 10^" eV; (ii) a local 
emissivity, which acts as an overall normalization 
factor on the UHECR spectrum; and (iii) a star for- 
mation rate (SFR). Here the underlying assump- 
tion is that the rate-density of GRBs follows star 
formation activity, expressed in terms of the mass 
processed into stars per comoving volume, and as- 
sumed to be traced primarily by hot young stars. 
Knowledge of the SFR is obtained by analyzing the 
blue/UV luminosity density using statistical sam- 
ples of galaxies. Blue light is thought to be a good 
proxy of star-formation activity, but extinction by 
dust complicates the measurement. 

Fig. 5 shows various SFRs based on different ap- 
proaches to the problem. SFR 4, from Hopkins and 
Beacom [82], is based on a compilation of IR, op- 
tical, and UV data. SFR 1 gives the lowest accept- 



able rate compatible with the data. SFR 3 gives 
the rate assuming large extinction corrections, and 
was used to fit the UHECR spectrum in [31]. With 
this rapidly increasing rate density of GRBs at red- 
shifts z w 1 - 3, a large pair production trough 
at w 10^^-^ eV is formed [83]. A rate density of 
GRBs following the extremely active SFR 5, based 
on IR luminosity density, is ruled out from calcu- 
lations of the UHECR spectrum, as can be seen in 
Fig. 6. 

A statistical study [84] of the redshift and 
opening-angle distribution of GRBs observed be- 
fore Swift, primarily Beppo-SAX GRBs, and the 
z-distribution of GRBs observed with Swift, re- 
jects SFR 3 for GRBs based on a standard rela- 
tivistic jet model. We concluded that the GRB rate 
density must monotonically increase to z w 5 - 7 
to explain the GRB statistics. The UHECR spec- 
trum from SFR 6, seen in Fig. 6, gives a reason- 
able fit to the HiRes data (the SFR 7 spectrum is 
nearly the same). Normalization to the Auger data 
would change the emissivity normalization down 
by a factor w 1/3. 

The idealized SFR oc (1 + z)^ sketched in Fig. 
5 was used by Berezinsky and his collaborators 
[85, 86] to describe the SFR activity of the UHECR 
sources, possibly blazars. Such extreme SFR ac- 
tivity produces a large pair production trough, and 
they proposed that the ankle/dip feature in the 
UHECR spectrum is due to photopair losses. Rea- 
sonable fits to the HiRes data were obtained with 
injection indices into intergalactic space s « 2.7. 



This model SFR oc (1 



can hardly be cor- 



rect, but models of blazars are more difficult than 
of GRBs by requiring both luminosity and density 
evolution. The long-duration GRB engine could 
very well be unchanged with epoch, but both due 
to black hole growth and fueling, the UHECR out- 
put from blazars must change with time. This be- 
havior is not well understood, so a mathematical 
model might be the best that can be done pending 
better studies (cf. [34]). 

Figure 7 shows predictions of the total diffuse neu- 
trino intensity spectrum for the various GRB mod- 
els in Figure 6, compared with the ANITA sensitiv- 
ity for a 45 day flight [87], and an estimated sen- 
sitivity for the proposed ARIANNA project [88]. 
Long-duration balloon-borne high-energy neutino 
telescopes can already start to test SFR 3 used in 
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Figure 7: Comparison of model calculations of 
GZK neutrino intensities for neutrinos of all fla- 
vors with sensitivity curves as labeled. 

our GRB model for UHECRs [31], and is close to 
discovering the guaranteed flux of GZK neutrinos. 
The cosmogenic neutrino flux from GRBs is diffi- 
cult to detect with IceCube, which is more sensi- 
tive to cosmic PeV sources (see [89] for IceCube 
sensitivities). 

Survival of UHECR Ions in the Blazar 
Environment 

Spurred by the Auger results on composition [5], 
we examine whether UHECRs can leave the blazar 
environment without being significantly photo- 
eroded. This assumes that blazar jets accelerate 
UHECRs, which is feasible given that the UHECR 
spectrum can be fit using simple assumptions for 
the SFR activity and injection spectrum [85, 86]. 7 
rays from blazars could be a consequence of ultra- 
relativistic hadrons in blazar jets [90, 68]. Cosmic 
rays accelerated in the inner blazar jet can power 
the knots and lobes of FR II radio galaxies through 
UHE beamed neutral production and escape due 
to photopion interactions of UHECR protons with 
ambient radiation in the inner jet [68]. The for- 
mation of a neutral beam in FR II radio galaxies 
could explain the Chandra X-ray emission from 
knots and hot spots in radio galaxies [91] by a sec- 
ond synchrotron component induced by UHECR 
activity [92]. 



Crucial to efficient neutrino and neutral beam 
power is that the blazar jet is found within the BLR 
where the scattered radiation field is intense. The 
question of the location of the 7-ray production site 
[93, 94], which would almost certainly be the loca- 
tion of the UHECR accelerator, will be settled by 
correlated GLAST/radio observations. 

Ions accelerated in the inner jets of BL Lac objects 
and FR I radio galaxies can pass through the broad 
line region (BLR) without significant photodisinte- 
gration, as we now show. Underlying this estimate 
is the unification scenario for blazars as expressed 
in [95], in which FR Is are the parent population of 
BL Lac objects, and FR lis are the parent popula- 
tion of FSRQs. 

Let the blazar BLR be approximated by a 
spherically-symmetric distribution of scattering 
electrons with density 7io{x) at distance x from the 
central source. The scattered radiation density in 
blazars can be estimated by noting that a fraction 
w no{x)aT{x/2) of the ambient photons will be 
scattered within a shell of width « x/2. For an 
isotropically emitting central source of radiation 
with photon production rate Npii{e) per unit di- 
mensionless photon energy e = hv/rrieC^, the am- 
bient photon density from the central source emis- 
sion is n,pfi{e;x) = Nph{e)/ATTx'^c. For assumed 
isotropic Thomson scattering, the spectral density 
of scattered radiation is 



=(e;a;) 



no{x)aTNph{e) 



(17) 



Sttxc 
[96, 97, 98]. 

The central source emission is assumed to be radi- 
ated by an accretion disk around the supermassive 
black hole. We represent the blue-bump emission, 
commonly observed in Seyfert galaxies (it is diffi- 
cult to observe in blazars because of the luminous 
jet radiation), by a Shakura-Sunyaev disk spectrum 
of the form 

e^^/'^ cxp(— e/e„ 



Nphie)^Lo 



(18) 



normalized to the total Shakura-Sunyaev disk lu- 
minosity Lq = 10'*^L46 ergs s~^, where r(4/3) = 
0.893, and emax is the maximum photon energy 
radiated in the disk spectrum. For the UV bump 
observed in supermassive black holes, e_5 = 
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Figure 8: Comparison of photo-nuclear destruction 
cross section for Fe from [12] and [101]. 



The blue bump in 3C 273 reaches a jyF,^ peak flux 
of « 3 X 10^^" ergs cm^^ s^^ at « 10 eV, corre- 
sponding to ;yL,^ w 2 X 10**^ ergs s~'^ and e_5 = 2 
[99]. The hard UV emission component observed 
from 3C 279 with lUE [100] has a peak vL^ lu- 
minosity « 3 X 10^'^ ergs s~^. Its effective tem- 
perature is w 20, 000 K, corresponding to a mean 
dimensionless photon energy « 10^^. 

A (5-function approximation for the photodisinte- 
gration cross section of a nucleus with atomic mass 

A = 56^56 is 



(19) 



[101], where ctq.a = mb, A ^ 15.6, 

6^,0 — 83. 5^°'^^, and is the invariant dimen- 
sionless photon energy in the ion's rest frame. A 
comparison of the photonuclear destruction cross 
sections for Fe is shown in Fig. 8. The (5-function 
approximation should be fairly good in all cases 
where the target photon spectrum is not too hard. 

For an ultra-relativistic ion passing through an 
isotropic radiation field, the probability per unit 
pathlength for the ion to photodisintegrate by in- 
teracting with ambient photons is given by 



(j)nuc 



dx 



(27')-' X 



roo /■27e 

/ de e~'^nph{e) / dcr Cr (TA{er) , (20) 
Jo Jo 

where the particle Lorentz factor 7 = E/AnipC^. 



The probability of an UHECR ion photo- 
disintegrating as it travels through the AGN BLR 
is, using eqs. (17) - (19) in eq. (20), simply 



P. 



(pnuc 



o-o,AA£r,oT^o(x)aT^o 
32r(4/3)72m3e3 



dx 



/ du u~^/^ exp(-«) , 

Jun 



max uq 

(21) 

where uo = erfl/2^emax- 

Here we take the typical extent of the BLR as « 
0.1 - 1 pc, and mean optical depth tt ^ 0.1. The 
BLR medium is probably clumped in rather dense 
clouds and has a strong density gradient from the 
BLR to the narrow Une region [102, 103], but here 
we approximate it as being rather uniform within 
a shell of radius lO^^i^ig cm with Thomson depth 
Tt = 10~2r_2, so that the mean BLR density is 
no = 1.5 X 104r_2/i?i8 cm'^. 
Eq. (21) is easily solved in the Umit uq <C 1 or 



4^ - eV . 

e-5 



The result is 



Pa. —012 
^ (pnuc — L/.j.^ 



pl/3 4/3p 
^20 ^-5 -^^18 



(22) 



forS20»10-3Al?6^7e-5. 
This result indicates that for typical parameters 
that may characterize the BLR of BL Lac objects, 
with T_2 ^ 1 and Lq w 10""^ ergs s^^, UHECR 
ions will escape without undergoing photonuclear 
breakup. The BLR environment may pose a haz- 
ard to UHECR ion escape in the luminous FR II 
radio galaxies and FSRQs with broad optical emis- 
sion lines. But note that our calculation only con- 
sidered a single interaction with the loss of one or 
a few nucleons from the UHECR ion. The corre- 
sponding probability for complete breakup will be 
a factor ^ A/2 smaller. 

This estimate suggests that UHECR ions can es- 
cape from BL Lacs and also from FSRQs, except 
in the cases of the most luminous blazars with thick 
columns of BLR material. The scattered radia- 
tion field in blazars is a convolution of the cen- 
tral source luminosity and surrounding gas distri- 
bution. The 77 attenuation process gives a sepa- 
rate probe of this radiation field [96, 104, 105]. By 
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Figure 9: Measurements of the EBL at optical 
and infrared frequencies [106], plotted in terms of 
spectral energy density eu(e), and fits to the EBL 
using modified blackbody functions. Data trian- 
gles pointed upwards refer to lower limits, and data 
triangles pointed downwards refer to upper limits. 
The spectral energy density of the CMBR at z = 
is also shown. 

jointly analyzing photodisintegration and 77 pro- 
cesses, GLAST data can be used to determine if 
the black hole jet environment limits UHECR es- 
cape. 

The Extragalactic Background Light 

Fig. 9 shows measurements of the intensity of the 
unresolved IR and optical EBL at the present epoch 
from the review by Hauser and Dwek [106, 107], 
including an upper limit at « 0.1 eV inferred from 
7-ray observations [108]. Motivated by the appear- 
ance of two distinct peaks in the SED of luminous 
infrared galaxies [109, 110], by synoptic spectra 
of galaxies in the Sloan Digital Sky Survey [111], 
by SEDs of nearby normal galaxies, including data 
from the Spitzer Space Telescope [112] and normal 
star-forming galaxies like the Milky Way [113], I 
fit these two peaks with modified blackbody func- 
tions. The lower energy emission feature peaking 
near w 0.01 eV, probably due to radiation repro- 



Table 2: Properties of the Dust and Two Stellar 
Components 



Component 


T(K) 


Uo 


k 






(IQ-^'^ ergs cm^^) 




Dust 


31 


0.273 


3.8 


Star 1 (HI EBL) 


7100 


1.1 


2.0 


Star 1 (LO EBL) 


7100 


1.1 2 


2.0 


Star 2 


16,600 


0.5 


3.0 



cessed by dust, is referred to as the dust compo- 
nent. The higher energy emission feature peaking 
near 2 eV is referred to as the stellar component. In 
our calculations, two modified blackbodies make 
up the stellar component. More terms can be added 
as required. 

The modified blackbody spectral energy density is 
written in the form 

eu{e) = Uo = meC^e^nphit) , 

exp(w) — 1 

(23) 

where w = e/8. For a blackbody, k — A and 
uq = 87rTOec2e^(l + zf/\l = 6. 37x 10-14(1 + 
zY ergs cm-'^. The fits to the data in Fig. 9 use 
the parameters given in Table 2. The high and low 
EBLs differ only by a factor 2 for the intensity of 
the lower temperature stellar radiation field. 

Fig. 10 shows the optical depth to 77 pair produc- 
tion attenuation for 7 rays with measured energies 
E detected from the TeV XBL lES 1101-232 at 
z = 0.186 [114, 115]. Separate components for the 
CMBR, dust, and stellar radiation fields are shown 
for the low EBL in the figure. In making this calcu- 
lation, only the CMBR field evolves with redshift, 
and the dust and stellar radiation field energy den- 
sities remain roughly constant. This assumption 
becomes increasingly less reliable at higher red- 
shifts. 

The attenuation factor, from which the intrinsic 
spectrum of lES 1101-232 is obtained, is plotted 
in Fig. 11. As can be seen from the index, the 
use of the low EBL means that the intrinsic pho- 
ton spectral index of lES 1101-232 from w 0.2 - 
3 TeV is « —2.0. If we adopt as a general rule, 
consistent with our knowledge of the GeV spec- 
tra of FSRQs [44] and GeV - TeV spectra of BL 
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Figure 10: Optical depth for a source at z — 0.186 
to 77 attenuation for the low and high forms for 
the EBL shown in Fig. 9. 

Lacs like Mrk 421 and Mrk 501, that the intrin- 
sic spectrum is softer than —2, then the low EBL 
is favored (cf. [116, 117]). A low EBL between 1 
and 10 ^ solves the problem of the unusually hard 
7-ray spectrum of lES 1101-232 [118, 114], and 
avoids having to construct acceleration scenarios 
not operating in Mrk 421 and Mrk 501, and to ex- 
plain lack of evidence of hard synchrotron compo- 
nents associated with a hard electron component in 
TeV/XBLs. 

The low EBL with a steep 2-10 micron spectrum 
favored to explain the TeV blazar data is in gen- 
eral agreement with the Primack model for galaxy 
formation [119], which considers star formation, 
supernova feedback and metal production in merg- 
ing dark matter halos. The low EBL disagrees 
with the EBL derived by Stecker and collaborators 
[120, 121, 122]. The problem is that their empiri- 
cal data base relies heavily on IRAS data at 12, 25, 
60, and 100 /i and uses a poor representation of the 
galactic SEDs between 1 and 10 /i. Individual IR 
and normal galaxies show far more structure in this 
region than considered. 




0.1 1 10 

E (TeV) 



Figure 1 1 : Attenuation factor for the low and high 
forms of the EBL, for lES 1101-232 at z = 0.186. 
Inset shows the effects of the low and high EBLs 
on the intrinsic spectrum of lES 1101-232 [118, 
114]. 

The 7-Ray Horizon 

The low and high EBL SEDs represent the likely 
range of the local z <C 1 IGM IR and optical radi- 
ation fields. We can use this field to calculate the 
photon horizon where t^^{E, z) — 1 [123, 124] in 
the limit z <C 1, when the IR and optical radiation 
fields have not changed appreciably over time. The 
result is shown in Fig. 12. This diagram is primar- 
ily illustrative, and in some respects misleading. At 
the redshift of lES 1101-232, namely z = 0.186, 
this diagram says that the exponential cutoff en- 
ergy is at « 300 - 400 GeV, and that the low and 
high EBLs are not significantly different. In fact, 
no exponential cutoff is seen in the lES 1101-232 
TeV spectrum (see inset to Figure 11), because the 
actual attenuation is very sensitive to the full spec- 
trum of the EBL. 

This figure does illustrate that PeV 7-rays can be 
detected from sources within our Galaxy, though 
they might be subject to modest attenuation from 
the CMBR. The IR and stellar radiation fields can 
also contribute significant 77 opacity at ~ 100 
TeV. Anisotropy effects of the radiation fields on 
opacity have recently been calculated by [113, 
125]. Attenuated spectra of TeV - PeV 7-ray 
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Figure 12: Low redshift {z <C 1) 7-ray horizon 
giving the relationship between photon energy and 
z or distance where r-y^ {E, z) = 1 for the low and 
high EBLs shown in Fig. 11. Calculated up to z = 
0.25. 



sources could in principle give a distance measure 
of specific sources in the Milky Way and nearby 
galaxies. 

The UHECR Ion GZK Radius 

We now use our low and high estimates of the EBL 
to calculate the ion GZK radius. First we show, in 
Figure 13, the energy loss and interaction mean- 
free paths for UHECR protons interacting with 
the combined EBL and CMBR, with components 
as shown. The dust component makes a minor, 
~ 10% contribution at « 8 x lO^'' eV, but could be 
somewhat larger for a modified fit to the FIR EBL. 
Other than that, the energy-loss mean free path is 
essentially given by the CMBR result (Figure 1). 
Note also that the ratio of the CMBR energy loss 
and scattering mean-free paths is w 5 on the low- 
energy wing where protons interact with the expo- 
nential Wien portion of the blackbody distribution. 
This ratio, arising from the 20% inelasticity for di- 
rect and resonance pion production, decreases at 
4 X 10^° eV due to the greater fraction of mul- 
tipion interactions at higher energies. 
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Figure 13: Energy loss mean- free paths for pro- 
tons in the combined EBL and CMBR radiation 
fields, including photopion, photopair, and expan- 
sion losses. The different low and high EBLs make 
little difference on the UHECR proton spectrum. 

Before calculating ion mean-free paths, it is worth 
mentioning how the energy-loss formula for pho- 
todisintegration is calculated. When a single nu- 
cleon is ejected, then an energy loss cx 1/Aof the 
original energy E happens (technically, provided 
that the nucleon ejection is isotropic in the nucleon 
frame), and for the ejection of two nucleons, an 
energy loss — 2E/A occurs. For multi-nucleon in- 
jection, an average factor is used, given by Puget 
et al. (1976) [12], = 3.6/ A for 10 < A < 22, and 
= A.349E/A for 23 < A < 56. A low-energy 
threshold of 10 MeV is used [13]. Obviously, after 
a single interaction, the original nucleonic identi- 
fication is changed, so that new sets of loss rates 
have to be used for the daughter particles. The 
calculated MFPs have only a generalized meaning, 
but provide inputs for accurate Monte Carlo simu- 
lations. 

Figure 14 gives various contributions of the CMBR 
and low and high EBLs to the effective energy-loss 
rate of UHECR Fe in an IGM radiation field at 
the present epoch. Also shown is the photodisin- 
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Figure 14: Photonuclear MFPs for UHECR Fe, in- 
cluding separate contributions to photodisintegra- 
tion, and total energy-loss and interaction rates. 

tegration interaction rate for the ejection of at least 
one nucleon. In effect, the energy-loss MFP gives 
the distance an Fe ion would have to travel to be 
broken up into mostly protons and neutrons. For 
E >6 X 10^9 eV, Fe only has to travel O(Mpc) 
before being transformed to lighter elements, and 
could hardly be seen in abundance in the UHE- 
CRs above this energy unless UHECR sources re- 
side in the neighborhood of our Galaxy, includ- 
ing M31 and our satellite galaxies (possible for a 
GRB origin of the UHECRs). Secondary nucleons 
with A ~ 56/2 would be more prevalent due to Fe 
photo-erosion. 

Figure 15 compares different assumptions for the 
photopion energy-loss rate on the total energy- 
loss rate of UHECR O in the CMBR and the low 
EBL. Without a detailed physical model, the pho- 
topion cross section should go oc A^^^ for a quasi- 
spherical nucleon, with a photopion energy-loss 
rate oc A^/"^, giving the maximum MFP. More re- 
alsitically, only one pion is produced with near 
threshold energy in the interaction, so the inelas- 
ticity would be cx A^^, and the photopion energy- 
loss rate cx A^^^"^, giving the minimum MFP. In 



Figure 15: Comparison of the effects on the total 
energy-loss MFP of UHECR O in the CMBR and 
low EBL field for different assumptions of the in- 
elasticity in inelastic photopion production. 

multi-pion production, the inelasticity should be 
larger than the minimum, so the "true" photopion 
energy loss rate should reside between these two 
extremes. The difference in either case is not sig- 
nificant below 10^^ eV, as can be seen. 

Proton and ion MFPs for energy loss in the com- 
bined CMBR and low EBL are shown in Figure 16. 
Most dramatic is the rapid decrease of the MFP be- 
tween 10^^ and 10^'' eV, precisely where the HiRes 
and Auger Observatories discover the spectral soft- 
ening in the UHECR energy spectrum [1,2]. Keep- 
ing in mind that the essential destruction of a nu- 
cleus of a given type is proportional to the scat- 
tering rate, and that the complete breakup propor- 
tional to the energy-loss rate, then this figure shows 
that structure within 100s of Mpc become visible 
above w 6 x 10^^ eV. At these energies, all ions 
will have undergone significant breakup, so that an 
original enhancement of Fe would be broken up to 
smaller Z. With increased statistics, higher energy 
events should be identified with even closer struc- 
tures. 
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Figure 16: Comparison of energy-loss MFPs of 
UHECR protons and ions in the combined low 
EBL and CMBR. 



Figure 16 reveals a number of interesting things. 
First, whereas UHECR protons have significant 
photopair losses to w 10^* eV, there is no change 
in ionic composition below w 5 x 10^® eV. At 
these energies, all ion losses are due to expansion. 
Above « 5x10^* eV, corresponding to the ankle or 
dip energy, photopair and photodestruction losses 
become roughly equally important. The dust com- 
ponent of the EBL (whose SED is not that accu- 
rately known) is important for the rapidly increas- 
ing energy-loss rate between w 5 x 10^^ - 6 x 10^^ 
eV. Above « 6 x 10^^ eV, photodisintegration by 
the CMBR starts to dominate, reducing the interac- 
tion length to tens of Mpc or less. UHECR Fe in- 
jected at /t lO'^" eV will, in short order, degrade to 
lighter nuclei at lower energies, which have larger 
cross sections at a given energy to degrade the nu- 
clei. 

By the accumulation of lighter Z material be- 
low the GZK energy, an ankle and cutoff may be 
formed through the injection of UHECR ions with 
large Z ^ 30- 56. The formation of this dip, seen 
in AGASA, HiRes and Auger data, must be ex- 



plained in any viable model for UHECRs. The ap- 
parently natural explanation of the ankle as a pair 
production trough resulting from UHECR proton 
injection [83, 31, 86, 85] might also find an expla- 
nation in GZK nucleonic physics. Detailed study 
requires a Monte Carlo simulation. 

IGM Magnetic Field 

We use the Auger results [15] to set lower limits on 
the mean IGM magnetic field B. The equation for 
the deflection angle from a source at distance d is 
[25, 126] 

d dZeB 



2 go 5-iid(100 Mpc) 

10 E2()\/Ninv 



(24) 



Here Ni^y = d/A ^ 1 is the number of inver- 
sions of the magnetic field, also expressed through 
the magnetic-field correlation length A. If the two 
UHECRs within 3° of Cen A were accelerated by 
the radio jets of Centaurus A, the measured 3° de- 
flection implies that 



20 



{E/Q X IQi^ eV) 



If the UHECRs originated in fact from the AGNs 
in the Veron-Cetty and Veron catalog [127], taking 

d = 75 Mpc and 9^ = 3° gives 



0.9 



[E/Q X 10^^ cV) 



Both of these values are reasonable, and possi- 
bly compatible depending on gradients in the IGM 
field. The IGM magnetic field energy density is 
a small fraction of the CMBR energy density or 
the EBL energy density. This shows the potential 
power of Auger for measuring IGM fields [128]. If 
these IGM fields are accurate, 7-ray pulse broad- 
ening could not be measured [129]. 

The equation for the time delay due to the deflec- 
tion of UHECRs from an impulsive source is 



35052^^^3(10 Mpc) .Z.2 
N~ ^10^ 

(25) 
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[25, 126]. Recurrent events over decades or 
shorter could be observed from UHECRs pow- 
ered by nearby GRBs at the distance of Cen A or 
GRB980425/SN1998bw (d ^ 36 Mpc) [130] in 
favorable circumstances. 

Galactic Cosmic Ray Astronomy is Un- 
likely 

It is great that Auger [15] has opened the field of 
extragalactic cosmic-ray astronomy, but the Pierre 
Auger Cosmic Ray Observatory (and its Northern 
Hemisphere counterpart, the Telescope Array, with 
half of Auger's effective area) may prove less use- 
ful for undertaking cosmic ray astronomy of galac- 
tic sources, simply because these sources are too 
weak. 

To detect clustering in the arrival directions of 
signal (the definition of an astronomy), then the 
source distance 

E20 



rf< \/NinvrL = 10\/iV~ 



(Z/10)B^G 



kpc ; 
(26) 

otherwise the particle trajectories would be hope- 
lessly scrambled. To accelerate particles with en- 
ergy E, equation (10), requires source luminosities 

Writing this expression as a limit on particle en- 
ergy gives, using eq. (26), the maximum source 
distance 



d « 20 



^Ninv (V1038 ergss-1) 



B 



PC (27) 



(Z drops out of the expression). 

In the Milky Way, where a large scale ordered field 
is measured [131], probably not a large number 
of inversions could occur over a distance 100 
pc, so I argue that Ninv is not large. In this 
case, there are not sufficiently luminous sources of 
nonthermal power close to the Solar system that 
could be observed by the cosmic rays accelerated 
at that source. For example, the 7-ray luminos- 
ity of Geminga, d ^ 140 pc, is w 3 x 10^^ 
ergs s-^ [132]. The shell-type SNR RX J1713.7- 
3946, if 1 kpc distant, releases « 10^^ ergs s~^ 



in nonthermal 7-ray energy and ^ lOx more in 
nonthermal synchrotron radiation [133|. The shell 
SNR RX 10852.0^622, with d ^ 200 pc, releases 
PS 3 X 10^2 ergs s^^ [134] between 1 and 10 TeV. 

This would explain the futile search for Galac- 
tic sources of cosmic rays [23, 135], which could 
only be detected from rare nearby SNe or hyper- 
energetic GRBs in our Galaxy [136, 137]. 

UHE Source Neutrinos from Blazars 
and GRBs 

Discussion of neutrino production from GRBs will 
be kept short because it was recently reviewed 
[138, 139]. Our basic approach is to use mea- 
sured flare nonthermal fluence to get apparent 
isotropic energy in the comoving frame within the 
uncertainty of the Doppler factor. For GRBs, we 
use observational constraints, including the burst 
rate and typical energy release, to normalize the 
mean (volume- and time-averaged) energy emis- 
sivity needed to power the UHECRs, from which 
the amount of energy that a typical GRB must re- 
lease in the form of nonthermal hadrons can be de- 
rived. Our results [31] showed that for an origin of 
UHECRs from long-duration GRBs, with an upper 
SFR 3 (Figure 5), long-duration GRB blast waves 
must be baryon-loaded by a factor fcR ^ 60 com- 
pared to the primary electron energy that is inferred 
from the X/7 GRB flux. 

The GRB neutrino fluences are to first order pro- 
portional to the electromagnetic X/7 radiation flu- 
ence from a GRB. Neutrino fluences expected in 
the coUapsar GRB scenario from a burst with pho- 
ton fluence ^rad = 3 X 10~*ergcm~^, were 
calculated [140] for 3 values of the Doppler fac- 
tor 6 from a GRB at redshift z = I (h = 65). 
For a GRB coUapsar-model calculation, we as- 
sumed that the prompt emission is contributed by 
^spk = 50 spikes with characteristic timescales 
tgpk — Is each, which defines the characteristic 
size (in the proper frame) of the emitting region 
associated with each individual spike through the 
relation R'^^^. ~ tspk5/{l + z). 

The numbers of muon neutrinos for IceCube pa- 
rameters that would be detected from a single 
coUapsar-type GRB with a baryon-loading factor 
fcR = 20 for 5 = 100, 200 and 300 are = 
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1.32, 0.105 and 0.016, respectively. For the large 
baryon load required for the proposed model of 
UHECRs, our calculations showed that 100 TeV - 
100 PeV neutrinos could be detected several times 
per year from all GRBs with km-scale neutrino de- 
tectors such as IceCube [140, 31]. Detection of 
even 1 or 2 neutrinos from GRBs with IceCube or a 
northern hemisphere neutrino detector will provide 
compelling support for a ORB origin of UHECRs. 
See [138] for more details. 

Detailed numerical simulations to calculate neu- 
trino production in colliding shell scenarios of 
GRBs, including the diffuse neutrino intensity, are 
given in [141]. Calculations of neutrino fluxes 
during X-ray flares found with Swift was treated 
in [142], and neutrinos and UHECRs from low- 
luminosity GRBs, in [143]. 

Following a similar methodology and normaUzing 
to blazar flare fluences, Armen Atoyan and I [9, 68] 
calculated neutrino fluxes from blazars; see [144] 
for a recent review. We found that for equal power 
into ultra-relativistic electrons, which makes the 
X/7 emission, and protons, IceCube could detect 
1 or 2 neutrinos from a powerful blazar flare Uke 
the 3C 279 flare in 1996 [78], of which GLAST 
should see one or two per month. 

Our principle result was that FSRQs, with strong 
scattered radiation field, provide much more ef- 
fective target photons for photopion production in 
their jets than do BL Lac objects. Thus powerful 
FSRQs like 3C 279 (z = 0.54), PKS 0528+134 
(z = 2.06), CTA 102 (z = 1.037), of 3C 454.3 
(z = 0.86) are a better targets for neutrino tele- 
scopes than nearby BL Lac objects. 

It will be important to check these conclusions with 
new, incoming 7-ray data. The incredibly short 
flaring timescale in the TeV XBL PKS 2155-305 at 
z = 0.116 on time scales of 100s of seconds [115], 
with bolometric apparent 7-ray powers of 
ergs s~^ (depending on EBL assumption), repre- 
sents a different regime for neutrino production in 
BL Lacs than considered earlier, and a probe of 
black holes on small size scales [145]. For blazar 
modeling of all types, GLAST will provide an ex- 
cellent data base to determine parameters used in 
models for neutrino production, (see [146, 147] 
for blazar neutrino searches.) 



Summary 

Auger has already contributed three major discov- 
eries to cosmic-ray physics: 

1 . the GZK cutoff, also found with HiRes; 

2. Mixed ionic composition in the UHECRs up 

to a few X 10^*^ eV; and 

3. statistical demonstration that the clustering 
of 27 UHECRs wifli energies > 6 x 10^^ 
eV follow the matter distribution as traced 
by nearby ( < 75 Mpc) AGNs. 

These results establish without doubt that UHE- 
CRs originate from astrophysical sources. With 
this information, and based on past theoretical 
work, I argue that GRBs and radio-loud AGNs, 
which are classified as blazars when viewed on- 
axis, are the most probable sources of UHECRs, 
as can be demonstrated by hadronic 7-ray signa- 
tures. The most convincing evidence would be 
direct detection of PeV neutrinos from UHECR 
sources with IceCube or KM3NET. Detection of 
GZK EeV neutrinos from photopion interactions 
of UHECRs as they propagate through EBL with 
an Askaryan telescope hke ANITA will also im- 
portantly test astrophysical models. 

In this paper, I used 7-ray observations to con- 
strain the spectrum of the EBL between 1 and 
100 /i, where it is poorly known. By adopting as 
a general rule, consistent with observations, that 
blazar 100 Me V- TeV 7-ray spectra are steeper 
than —2, a low EBL, shown in Figure 9, was fa- 
vored. With this EBL, we can deconvolve the 
intrinsic source spectrum of low-redshift sources, 
and calculate the photodisintegration rate of differ- 
ent UHECR ions. The GZK curves for different 
ionic species obtained with the low EBL is shown 
in Figure 16. These results are in accord with the 
location of the GZK spectral cutoff measured with 
Auger and HiRes, and a mixed ionic composition, 
because protons are more difficult to accelerate, 
and Fe wiU be broken up. 

GRBs and blazars are both viable sites for UHECR 
acceleration, consistent with the Auger results. 
What will be important is to associate UHECR ar- 
rival directions with a subset of AGNs or galax- 
ies. Small metal-poor bluish star-forming galax- 
ies are likely hosts of long-duration GRBs [148]. 
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Radio galaxies would host misaligned blazars, and 
slightly misaligned blazars could be 7-ray dim and 
cosmic-ray bright. GLAST will be able to see 7- 
ray dim AGNs (the radio galaxies M87 [149], Cen 
A [54J, and NGC 6251 [150J are definite, likely, 
and probable 7-ray sources, respectively). Asso- 
ciating radio galaxies or GLAST 7-ray galaxies in 
the nearby universe with arrival directions of UHE- 
CRs could decide the question of blazar origin of 
the UHECRs. With the 7-ray, cosmic ray, and neu- 
trino observations, it is likely that the problem of 
UHECR origin will soon be solved''. 
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